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Abstract 

In this paper we consider a branching particle system consisting of particles moving according to 
the Ornstein-Uhlenbeck process in R'' and undergoing a binary, supercritical branching with a constant 
rate A > 0. This system is known to fulfil a law of large numbers (under exponential scaling). In the 
paper we prove the corresponding central limit theorem. The limit and the CLT normalisation fall 
into three qualitatively different classes. In, what we call, the small branching rate case the situation 
resembles the classical one. The weak limit is Gaussian and normalisation is the square root of the 
size of the system. In the critical case the limit is still Gaussian, however the normalisation requires 
an additional term. Finally, when branching has large rate the situation is completely different. The 
limit is no longer Gaussian, the normalisation is substantially larger than the classical one and the 
convergence holds in probability. 

We prove also that the spatial fluctuations are asymptotically independent of the fluctuations of the 
total number of particles (which is a Galton- Watson process). 

MSG: primary 60F05; 60J80 secondary 60G20 
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1 Introduction 

We consider a branching particle system {Xt}^^^ as follows. The system starts off at time t = from a 
single particle located at x G M*^. The particle moves according to the Ornstein-Uhlenbeck process in R'^ 
and branches after exponential time with parameter A > 0. The branching is supercritical and given by 
the generating function 

F{s) ■.= ps^ + {l-p), p> -. 
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The offspring particles follow the same dynamics. This system will be referred to as the OU branching 
process. Formally, the system is identified with the empirical process, i.e. X is a measure-valued process 
such that for a Borel set A, Xt{A) is the (random) number of particles at time t in A. We recall that the 
Ornstein-Uhlenbeck process is a time homogenous Markov process with the infinitesimal operator 

L := ^cr^A-^xoV. (1) 

where cr, /i > and o denotes the standard scalar product. The Ornstein-Uhlenbeck process has a unique 
equilibrium measure ip (to be described later). 

Systems of this type may be regarded as consisting of two components, namely the genealogy part and 
diffusion part (for this reason they are sometimes called "branching diffusions"). The genealogy part, being 
the celebrated Galton- Watson process is well-studied. In our paper, the expected number of the progeny 
of a particle is strictly greater than 1, therefore the system is supercritical. The expected total number of 
particles grows exponentially at the rate 

\p := {2p - 1)A. (2) 

After a long time the positions of two "randomly picked" particles are "almost independent" random 
variables, which suggests the following law of large numbers 

^ ((^,/) W, a.s. (3) 

where Exf^ is the event that the system does not become extinct, \Xt\ denotes the number of particles at 
time t and / is a bounded, continuous function. This is indeed the case as follows from |17| . In Theorem 
13.11 we obtain ([3]) for a slightly more general class of functions. This is, however, only a preparatory 
step for our main goal which is the corresponding central limit theorem. The second order behaviour 
depends qualitatively on the sign of Ap — 2/i. Roughly speaking this condition reflects the interplay of two 
antagonistic forces, the growth which is local and makes the system more coarse and the movement which 
tends to smooth the system (higher /i implies "a stronger attraction of particles towards 0"). Now we are 
going to describe the behaviour of the spatial fluctuations: 

F,-\{XtJ)-\Xt\{ipJ)), (4) 

where Ft is some not necessarily deterministic norming. We will describe the situation on the set of 
non-extinction Ext^. 

Small branching rate: Ap < 2/i. Our main result is contained in Theorem 13.21 In this the case "the 
movement part prevails" and the result resembles the standard CLT. The normalisation is given by 
Ft = \Xt\^^'^ (which is of order e^'*'''/^^*). Moreover, we obtain the limit which is Gaussian (though 
its variance is given by a complicated formula) and does not depend on the starting position x. Let 
us also note that a random normalisation is quite natural as it "filters" out the fluctuation of the 
total number of particles. 
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Critical branching rate: Ap = 2/i. Our main result is contained in Theorem 13.61 In this case "the 
branching prevails". The behaviour of the fluctuations slightly diverge from the classical setting. 
The normalisation is bigger: Ft = t^^'^\Xt\^^'^ . The limit still does not depend on the starting 
condition and is Gaussian but its variance depends on the derivatives of /. To explain this we notice 
that the branching is so fast that the fluctuations are are not smoothen by the motion and become 
essentially local. In consequence they give rise to a spatial white noise and larger normalisation is 
required. 

Large branching rate: Ap > 2^. Our main result is contained in Theorem 13.101 In this case not only 
does the branching "prevail" but also "the motion badly fails to make any smoothing". The normal- 
isation is even bigger: Ft = e^^^~^^^ and we have Ap — /x > Ap/2. The limit is no longer Gaussian, it 
is given by (/, grad ip) o J (where J is the limit of a certain martingale). What is perhaps surprising, 
the limit holds in probability. The first term, (/, grad(/3), means that alike the critical situation the 
branching is fast enough to produce some sort of a white noise. Even more, it is so fast that the 
limit depends on the starting condition and in fact, up to some extent, the system "remembers its 
whole evolution", which is encoded in J. 

In either case we prove also that the spatial fluctuations become independent of fluctuations of the total 
number of particles as time increases. 

To our best knowledge so far there have been no CLT-type results of a similar flavour in the fleld of 
branching diffusions and the research effort was concentrated instead on proving laws of large numbers 
for more and more general branching systems. The only CLT result we are aware of is contained in [SJ 
Proposition 6.4]. Their setting is somewhat different as they consider "running fluctuations" of the form 
F^r^ {{Xt+T, f) — {Xt, ft,T)) as T — )• +00, where is a normalisation and fT^t is a certain transformation 
of /. Besides many advantages their approach fails to capture the emergence of three qualitatively different 
cases described above. More detailed discussion is contained in Remark 13.151 It is also noteworthy that 
our results open a possibility of further research. We list a few most promising possibilities in Remark 
13.161 and Remark 13.171 

At this moment we would like to announce our parallel paper [1]. In [1] we consider the [/-statistics 
of the OU branching system, namely expressions of the form 

Unf)-= E fiXtiii),Xtii2),...,Xtiin)). (5) 

n,j2,.--*n = l 

ik¥=ij if k^j 

We obtain both the law of large numbers (i.e. an analogue of (|3|)) and CLTs, which also fall into three 
categories corresponding to the cases described above. Moreover, at this point we also advertise a forth- 
coming work of the second named author |19) , which is devoted to studies of the CLT for superprocesses 
based on the Ornstein-Uhlenbeck process. Qualitatively the results of |19| are the the same as the ones 
presented in this paper. 
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Our proofs utilise a mixture of techniques used for the branching particle systems (e.g. the Laplace 
transform and the log-Laplace equation, coupling and decoupling). Although the proof schemes loosely 
resemble known techniques (e.g. are similar to proofs in [U Section 1.13]) they required many improve- 
ments. In our proofs we used also the fact that the Ornstein-Uhlenbeck process has a particularly explicit 
and traceable structure. Improving this part seems to be an interesting research problem. 
The studies of branching models of various types have a long history, we refer the reader to [H [TTl HI [7] (the 
list is by no means exhaustive). It has been known for a long time that the behaviour of branching systems 
differs qualitatively for (sub)critical and supercritical cases. The former become extinct almost surely, their 
limit properties are studied after conditioning on non-extinction event (e.g. Yoglom type theorems) or as 
a part of larger infinite structures (e.g. Galton- Watson forests, random snakes, continuum trees e.g. |13|). 
The latter grow exponentially fast (on the set of non-extinction), which makes them possible to be studied 
using laws of large numbers, starting with the celebrated Kesten-Stigum theorem (|U Theorem 1. 10.1]). 
Such theorems were also proved for branching particle systems, they go back to [2l [3] and more recently 
[lOj (which was the main inspiration for our paper). It presents a law of large number for a large class 
of supercritical branching diffusions and admits unbounded space-dependent branching intensity. Also we 
would like to mention again paper [5] in which systems with non-local branching (i.e. particles may jump 
upon an event of branching) were studied. The article [S] presents a law of large numbers as well as a 
central limit theorem, though in a different spirit than ours (see comparison in Remark I3.15p . Due to 
its excessive introduction |5] is also an excellent resource of biological motivations for study of branching 
diffusions. 

The article is organised as follows. The next section presents notation and basic facts required further. 
Section [3] is devoted to presentation of results. Proofs are deferred to Section 3] and the Appendix. 

2 Definitions and notation 
2.1 Notation 

For a branching system {X^j^^Q, we denote by \Xt\ the number of particles at time t, and by Xt{i) - the 
position of the i-th. (in a certain ordering) particle at time t. Typically, we use or P^^ to stress the 
fact that we calculate the expectation for the system starting from a particle located at x. Sometimes we 
use also E and P when this location is not relevant (e.g. if we calculate the number of particles in the 
system). We will refer to system starting from a single particle at time t = located at x shortly, as 
the OU branching system starting from x G M'^. 
For a function / : — )• R, we will denote 

{XtJ)=Y,fiM^))■ 

1=1 
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By —7- we denote the convergence in law. We use to denote the situation when an equahty or 

inequahty holds with a constant c > 0, which is irrelevant for calculations. E.g. /(x) ~ g{x) means that 
there exists a constant c > such that f{x) = cg{x). 

Let X o y = Yli=i^iyi denote the standard scalar product oi x,y G W^. Moreover, = -^x o x is the 
standard Euclidean norm in M'^. 

We use also {f,fi) := J^^ f{x)fj.{dx). 

In the paper we will use the space 

V = ViW^) := |/ : M"* H> M : / is continuous and 3.„ such that as ||x|| +oo| , 

that is the space of continuous functions which grow at most polynomially. 

Given a function / G V{W^) we will implicitly understand its derivatives (e.g. ^) in the space of the 
tempered distribution (see e.g. [20| p. 173]). 

Basic facts on the Galton- Watson process The number of particles {|Xj|}^~^q is the celebrated 
Galton- Watson process. We present basic properties of this process used in the paper. The main reference 
in this section is [3]. We already introduced the growth rate ([2|) (e.g. [H Section 1.6]). The process 
becomes extinct with the probability pe := (see [H Theorem 1.5.1]). We will denote the extinction 
and non-extinction events by Ext and Exf^ respectively. The process 

Vt := e~^-'\Xt\ 

is a positive martingale. Therefore it converges (see also [U Theorem 1.6.1]) 

— ^ ^oo) CL-S- as i — )■ +00. (6) 

Fact 2.1. We have {Kx> = 0} = Ext and conditioned on non- extinction Voo has the exponential distribution 
with parameter . We have ¥,{Voo) = 1 and Var(yoo) = 2p-i ■ IEe~^'^p*|Xf |^ is uniformly bounded, i.e. 
there exists C > such that for any t > we have Kr^e~'^^p^\Xi\'^ < C. Moreover, all moments are finite, 
i.e. for any n G N and t > we have E|Xf|'" < +00. 

The proof is deferred to the Appendix. We will denote the variable Voo conditioned on non-extinction 
by W. The density of the invariant measure of the Ornstein-Uhlenbeck process defined by ([1]) is given by 

3 Results 

This section is devoted to the presentation of our results. The proofs are deferred to Section |4| Our first 
aim is to present a central limit theorem corresponding to the following law of large numbers (closely 
related to [101 Theorem 6] or [5l Theorem 4.2]) 
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Theorem 3.1. Let {Xt}^^Q be the OU branching system starting from x G M.'^. Let us assume that 
f G V{R'^). Then 

lim e""^*"* {Xt, f) = (/, (f) Voo in probability, 

t— S- + CXD 

or equivalently on the set of non- extinction, Ext^, we have 



lim \Xt\ ^ {Xt, f) = {f, If) in probability. (8) 

t— >-+oo 



Moreover, if f is bounded then the almost sure convergence holds. 

3.1 Small branching rate: Ap < 2/i 

We denote f{x) := f{x) — {f,<f) and 



a 



f 



^, p) +2Xp [ ^ (v, (e^'-/'>Tsf)) ds. (9) 







The above expression becomes perhaps less cryptic when expressed in the base of the Hermite polynomi- 
als. Let {i^i},j>o be the probabilistic Hermite polynomials (see. e.g. |18| page 5]) and let fii,i2,...,i^ '■= 
f^d f{x)Ylj^-i^Hi^{xj)ip{x)dx. Obviously, /o,o,...,o = 0. Moreover, we have 



a 



f 



We will comment on this formula in Remark 13.141 

Let us also recall ^ and that W is Voo conditioned on Ext^. The main result of this section is 

Theorem 3.2. Let {Xj}^>q be the OU branching system starting from x G M'^. Let us assume Ap < 2^ and 
f G V{W^). Then a'j < +oo and conditionally on the set of non- extinction Exf^ there is the convergence 

where Gi ~ M{0, l/{2p — 1)), G2 A/'(0, a'j) and W, Gi, G2 are independent random variables. 

Remark 3.3. As we already mentioned in the Introduction this is the most classical case. One can check 
that the theorem is still valid when the first particle is distributed according to ip. The third term in the 
theorem can be written as: Sl=i ~ (/i 9^))' moreover (by our special assumption about 

the starting condition) each Xt{i) ~ (p. If these random variables were independent then the above sum 
would converge to AA(0, fjj), where a'j = /^^ (the random number of elements in the sum is only a 
minor obstacle). Therefore, the additional integral term in ([9]) reflects the dependence between Xt(f)'s. 

Remark 3.4. An important feature of our result is the factorisation of the fluctuations of the total mass 
process {|Xt|}j and the spatial fluctuations process i.e. {{Xt, f)}^. Using this fact we can easily prove a 
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central limit theorem corresponding to ([8]) with deterministic normalisation. On the set of non-extinction 
Exf^ we have 

e-V/2 ((Xi,/) - \Xt\ if, if)) ^'^ G2VW, 
where W, G2 are the same as in Theorem 13.21 

Remark 3.5. The convergence of the spatial fluctuations can also be regarded as convergence of random 
fields. It is technically convenient to embed the space of point measures into the space of tempered 
distributions 5'(M'^) (i.e. the dual of the space of rapidly decreasing functions 5(M'^)). The Levy continuity 
theorem holds in nuclear spaces (e.g. |14| Theorem B.]), hence the following 5'(R'^)-valued random variable: 

._ Xt - \Xt\ip{x)dx 

^* ^^v^ ' ^^^^ 

converges to a Gaussian random field M with covariance structure given by 

Cov((M, /i) , (M, h)) = {if, hh) + 2Ap (e(^-/2)-r./i) (e^^^'^^Tsh) ) ds, 

J 

where ]i(x) = fi{x) - (ifji). 

3.2 Critical branching rate: Ap = 2/x 

We denote 



ambiguity. Analogously as in (jlOp we can expresses a'j nicely using the Hermite expansion: 



2 4Ap^ 



i=l 

Note that the same symbol has already been used to denote the asymptotic variance in the small 
branching case. However, since these cases will always be treated separately, this should not lead to 

= (/i,o,...,o + /o,i,...,o + /o,o,-,i) ■ (12) 
Let us also recall ([6]) and that W is Voo conditioned on Ext^. The main result of this section is 

Theorem 3.6. Let {Xt}^^^ be the OU branching system starting from x G W^. Let us assume that 
\p = 2fi and f € V{W^). Then < +00 and conditionally on the set of non- extinction Exf^ there is the 
convergence 



where Gi ~ A/'(0, l/{2p — 1)), G2 ~ A/'(0, CTj) and W, Gi, G2 are independent random variables. 

Remark 3.7. We continue discussion from Remark 13.31 This time the theorem is less classical as the 
normalisation is larger. We interpret this as the fact that Xj(i)'s become more dependent as Ap increases 
relatively to ^. In other words the branching is so fast that any particle has many relatives which are still 
close to it. This also explains, at least on an intuitive level, the appearance of derivative in (|12p . One can 
think that the limit, in a sense of random fields, is a Gaussian white noise. 
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3.3 Large branching rate: Xp > 2fi 

We introduce the process 

\Xt\ 

Ht := e(-^''+^)* ^Xt{i), t> 0. (14) 

i=l 

We have 

Fact 3.8. Let Xp > 2/i. Then Ht is a martingale with respect to the filtration of the OU branching system. 
We have sup^ ^xHf < +oo, therefore there exists 

Hoo ■■= hm Ht, 

t— >H-oo 

where the convergence holds a.s. and in L^. When the OU branching system starts from 0, then the 
martingales Vt and Ht are orthogonal. 

The distribution of H^o depends on the starting conditions. 

Fact 3.9. Let {Xt}t>Q o,nd {Xt}t>o be two OU branching systems, the first one starting from and the 
second one from x. Let us denote the limit of the corresponding martingales by H^o, H^o respectively. Then 

where Voo is given by ([6]) for the system X . 

Let us denote by J the random variable i^oo conditioned on Ext^. 

Theorem 3.10. Let {Xt}tyo branching system starting from x G W^. Let us assume that 

Xp > 2/z and f G V{M.'^). Then conditionally on the set of non- extinction Exf^ there is the convergence 

fe-V|x,|. l^'|-^;'^- . <^"^>-f-l (,.F,G.(grad/.,>o J). (15) 
\^ y/\Xt\ exp((A-/i)t) j 

where G ^ M{0, l/(2p — 1)) and (W, J), G are independent. Moreover 

^~>^vt\Xtl ^^tJ)-\Xt\{f,^) \ ^ (grad/,v9) o H^), m probability. (16) 
exp((A-//)t) J 

Remark 3.11. As we noted in the Introduction this case hardly resembles the classical CLT. The con- 
vergence in probability is perhaps its most unexpected feature. Unfortunately, we cannot present any 
satisfying intuitive explanation. This phenomenon seems to be closely related to the fact that the branch- 
ing is so fast that the system is "not able to forget" the starting condition and in fact, up to some degree, 
it "remembers" its whole evolution (encoded in martingale Ht). Analogously to the critical case the limit, 
treated as a random field, it is some sort of white noise (not Gaussian). 
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Remark 3.12. We were not able to derive any explicit formula for the law of Hod- However we calculated 
some of its moments. As the formulas become lengthy, we assume that p = 1 and present only: 



^2 _ o T,. . IT ^4 _ 967' (16 + 397 + 307^ + 87=^) 



EoiH^y = 27, EoiH^Y 



9 + 277 + 2672 + 873 

14407^ (^36847 + 2856757 + 9480127^ + 17604207^ + 20054087*' + 14411207^ + 6421127"^ + 1635847'' + 184327*) 



(1 + 7)^(3 + 27){5 + 47)(5 + 67)(5 + 87) (6 + I77 + 1272) 

where 7 = (Ap//U — 2)^-*^. The even moments are as the distribution is symmetric. One can now check 
that Hoc is never Gaussian. Moreover, V^o and H^o are not independent (even though uncorrelated!). 
Their dependence is not trivial. By similar calculations of the fourth moment of xVoo + H^q we also 
checked that H^o is not of the form VoqG, where G is some random variable (not necessarily normal) 
independent of Voo- 

Remark 3.13. We suspect that the convergence in (|16|) is in fact almost sure. 
General remarks 

Now we will present general remarks common for all cases. 

Remark 3.14. The forms of (jlOp and (jl3p as well as the one of (which is nothing else than e~^^p~^^* {Xt,Hi)) 
seem to stem from some hidden underlying structure. This indeed is the case which, for notation reasons, 
we will explain for d = 1. Formally, one can write 

+00 

1=0 

where fi is the expansion of / in the Hermite base and for i >0 processes {-^ij^^o given by 

Mi = e-(^p-^^)* {Xt,Hi). 

One easily checks that M* are martingales. Understanding this expansion and its relation to our results 
will be subject of further research. An analogous expansion for a slightly different model was presented in 

m- 

Remark 3.15. A result close to ours was presented in [H Section 5.1], where the authors consider a model 
{y^l^^Q in which particles move according to a diffusion on the real line with the generator Lf{x) = 
b{x)f'{x) + ^—^f"{x). Upon the event of branching the position of mother (interpreted in [5] as its 
value) is split between its progeny according to a certain random law. Further, they define, a family of 
measure- valued stochastic process {'nT}t>(j 

/„T .\ ._ (Ap(i+r))/2 ( jXt+Tj) _ {XT,Ttf) \ , . 

where Ttf = e'^p^E^ {f,Xt). In a carefully chosen functional space (we skip the details for the sake of 
brevity) rj'^ — t-"^ t], where i] is a solution of the following stochastic equation 

{riu f) = f f {Lf{x) + J/(x)) r?,(dx)ds + VwWt{f), 
Jo 
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where W is an analogue of our Voo, W is a certain Gaussian martingale with values in a functional 
space and J is a certain operator related to the branching. As we already mentioned in the Introduction 
the expression pT|) represents "running fluctuations" i.e. study the differences between two time points 
at a finite distance t while our approach is more like studying difference between time t and infinity. 
Consequently, it does not seem to be a direct relation between our result and that of [5j. This is also the 
reason why their result does not change when the branching intensity is being altered. 

Remark 3.16. The most important extension of the present work will be to study systems with particles 
movement given by more general Markov processes (starting with diffusions). It is also very interesting to 
study the case of systems with non- homogenous branching rate (like the ones in |10|). These questions are 
harder to answer than our present result. We expect that spectral theory of operators will play a crucial 
role. One may suspect that formulation of the corresponding results will be much more involved comparing 
to the Ornstein-Uhlenbeck case. The exponential rate of convergence to the equilibrium measure is essential 
for a behaviour similar to our slow and critical branching case. This entails that the strange behaviour of 
the supercritical case may be common for many natural examples, e.g. the branching Brownian motion. 
In this paper an essential idea was to use an explicit coupling of Ornstein-Uhlenbeck processes. This 
methods has an advantage of conceptualising and simplifying the proofs. However, it is not amenable to 
generalisations. The first step to overcome this problem was made in |19| . Working with a similar model, 
the author was able to prove CLT results using only analytical properties of the Ornstein-Uhlenbeck 
semigroup. We believe that such approach will be much easier to extend. We expect that the task will be 
easiest in the "small branching rate" case. 

To give the reader a glimpse of forthcoming complications let us consider the branching particle system 
X in M? such that particles move according to independent Ornstein-Uhlenbeck processes with different 
drift parameters < ^2- Let us denote by ipi, ip2 the invariant measures on the first and second coordinate 
respectively. We fix some g\^gi : M i— )■ M such that {gi,ipi) = and {g2,'^2) = 0. Further we denote 
/i '■= gi ^ ^, /2 = 1 ® 52- We now consider 

(XtJi) = (X^gi) , {Xtj2) = {Xlg2) . 

In the above we stick to notation that X^{i),X^{i) denotes the first, respectively second coordinate of 
the position of the i-th particle at time t. The situation is "standard" when Xp < 2fii = 2 min(/ii, /i2)- 
By Theorem 13.21 in both cases we obtain convergence to a Gaussian limit (with normalisation 
This is no longer true when Xp > min(/ii, using Theorem 13.61 and 13.101 one easily checks that the 
normalisation "required" by {Xt, fi) is strictly larger than the one corresponding to {Xt, /2). Even in this 
simple example the situation becomes "singular". 

We expect that the phenomena described above will hold in more general cases. Therefore, in the 
subcritical case it should be relatively easy to obtain analogs of our results for more general motion 
processes. Contrarily, the supercritical case may turn out to be much harder and the description of the 
limit may depend on the motion (or its infinitesimal operator) in a complex way. 
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Remark 3.17. Now we list other possible extensions of our results. The first has already been obtained, a 
parallel paper [1] contains the corresponding results for [/-statistics (as it was described in Introduction). 
Secondly, corresponding results was also obtained for superprocesses in |19) . Another natural extension 
would be to consider more general branching laws. The first problem to consider is a law of large numbers. 
The statement of the Kesten-Stigum theorem suggest the weakest possible conditions and the Seneta- 
Heyde theorem suggests further extensions (we refer the reader to [6] and references therein). Secondly, 
we suspect that at least some of our CLT results hold when the branching law has finite variance. The 
case of infinite variance branching laws is interesting as well, one should expect a stable law in the limit. 
Finally, other interesting lines of research are the large deviation principle and functional convergence. 

4 Proofs 

4.1 Ornstein-Uhlenbeck process 

The semigroup of the Ornstein-Uhlenbeck process will be denoted by T- It can calculated with the following 
formula 

rtf{x) = {gt*f){xt), xt:=e~'''x, (18) 

where 

We denote ou{t) := Vl — e"^/^* and let G ^ (p. The semigroup Thas the following useful representations 

Ttfix) = [ f{xt - y)gt{y)dy = [ f (xe"^* + ou{t)y) ^{y)dy = E/(xe-'^* + ou{t)G). (19) 

The following coupling will be very useful for further analysis. 

Fact 4.1. There exists a probability space and two Ornstein-Uhlenbeck processes {i]t}t>o, {7t}t>o defined 
on this space such that r]Q = x and 70 = which fulfil 

111-11= xe~^\ a.s. 

Proof. Let rj be the (unique, strong) solution of the following stochastic differential equation 

d?7i = crd/3t - /xr/tdt, % = x, 
where /? is the standard Wiener process. We construct 7 by using the same (3 i.e. 

d7f = adj3t - n^tdt, 70 = 0. 

The result follows by the subtraction 

d(?7t - 7t) = -Kvt - lt)dt, r^o - 70 = x. 

□ 
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4.2 Laplace transform and moments 

Now we compute the Laplace transform of X and derive moments. We use standard techniques in the 
branching processes theory (see e.g. |15| ) hence we skip some details. Let / : M'^ i— )■ M+ be a continuous 
compactly supported function. We denote 

w{x,t,e) :=E,exp(-(0/,Xt)). 

By standard conditioning, renewal arguments and the branching property we obtain 

w{x, t, 9) = e-^'Tte~^f (x) + \ f e~^^'~''>Tt^sF{w{; s, 9)){x)ds, 



where F is the generating function of the branching law. The last expression writes as 
d 



dt 



w{x, t, 9) = {L- \)w{x, t, 6) + {Xp)w{x, t, ey + A(l - p), w{x, 0, i 



Hence the Laplace transform satisfies 

w{x, t, 9) = Tte-^^{x) + X [ Tt-s [pw\; s, 9) - w{; s, 9) + {l- p)] {x)ds. 

Jo 

We differentiate (j2ip with respect to 9. For A: > 1 we have 
w('^Hx,t,9) ={-^fTt \f\-)e-'f{-)\ (x) 



(20) 



(21) 



t 

+ A / Tt-s 
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<^(''\w^^\.,s,9)w^''-^\.,s,9)-w^'^\.,s, 

1=0 ^ ^ 



(x)ds. 



Note that this differentiation is valid by Fact 12. II and properties of the Laplace transform (e.g. |121 Chapter 
XIIL2]). We evaluate this expression at 9 = 



w^^\x,t,Q) ={-lfTtf\x) 



+ X Tt- 

10 



> J2 (^) ^^'^ (•> 0)w^''-'^ i; s, 0) + i2p - l)u;W (., s, 0) 



(x)ds. 



We recall that Xp = (2p — 1)A. One easily checks that 



^0 



.|;Q^.«(.,.,0)u;('=-')(.,.,0) 



By the properties of the Laplace transform the moments are given by 

E,((/,Xi))" = (-l)V")(x,i,0). 



(x)ds. (22) 



(23) 



By the formula above w^'^\x, t, 0) can be made meaningful also in situations when the Laplace transform 
is not well-defined. In particular using standard techniques one can show that (j22p is valid for f V (in 
fact in this paper we will need only moments up to order 4). 
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4.3 Weak convergence facts 

In this section we gather simple facts concerning the weak convergence. Let us denote by || • \\rpy the total 
variation norm on the set of probability measures. We have a simple lemma 

Lemma 4.2. Let P) he a probability space and Ai <Z A2 ^ F be such that ¥ {Ai) > 0. Let X be a 

random variable and ui, vi be its law conditioned on Ai,A2 respectively. Then 

P(^2)-P(Ai) 

Proof. Let i? be a Borel set, then 

'{{XeB}nAi) F{{XeB}nA2) 



'(X £ BlAi) -¥(X e B\A2) , 

' {{X eB}n Ai) (P {A2) - P (^1)) + (P {{X eB}nAi)-F {{x eB}n A2))F [Ai] 



<2- 



P(^l)P(A2) 

'{A2)-F{Ai) 



IP(^i) 

The conclusion holds by the fact that B is arbitrary. □ 

Let /-ii,/U2 be two probability measures on M, and Lip(l) be the space of continuous functions M 1— >• 
[— 1, 1] with the Lipschitz constant smaller or equal to 1. We define 

m{^ii,^X2)-= sup - (5,^2) I- (24) 

gGLip(l) 

It is well known that m is a distance equivalent to weak convergence (see e.g. [HI Theorem 11.3.3]). One 
easily checks that when /Ui,/X2 correspond to two random variables Xi,X2 on the same probability space 
then we have 



m{^ll,^Ji2) < ||Xi -X2II1 < ^J\\Xl-X2\\2. (25) 
4.4 Rate of convergence to invariant measure and approximations 

We will need also estimations of the speed of convergence to the invariant measure. Throughout the proofs 
we will denote 

f-=f-{f,v>)- 

In proofs below it will be convenient to have some additional regularity conditions. Given f & V and 
14 > we denote 

Ux) = TuRx). (26) 

One easily checks that Z„ is a C°° function and its derivatives grow at most polynomially. Moreover, we 
have {lu, '■p) = 0. 



13 



Lemma 4.3. Let f G V{M.'^) and f be defined as above. Then there exist constants C,n > such that for 
any t > 1 we have 

Ttf{x) < C(l + IkDe-'^*, 71/(0) < Ce-^/^*. (27) 



Moreover 



lim e^*Ttf{x) = xo (grad /, (p) , 



where grad / is understood in a weak sense. There exist C, h such that 

e'^^TJix) -xo (grad /, if) < C{1 + \\xf)e~''K (28) 
Moreover, there exists a function c : M+ i— >■ M+ such that c{u) — )■ as u \ and 

Ttilu - f){x) < c(n)(l + llxDe-'^*. (29) 
Proof. Let us recall equation (jl9p . We write 

e^(*+")7^+n/(x) = e^"e^*7I/n(x) = e'^" / e'^*(Z„(xe-'^* + ou{t)y) - Uy)My)dy. (30) 



Using the mean value theorem we get 

h{x, y, t) := e^*(/„(xe-^* + ou{t)y) - l^{y)) = {x + e^'\ou{t) - l)y) o grad /^(xo), (31) 

where xq is some point on the interval joining y and xe~^* + ou{t)y. From this representation we obtain 
\h{x,y,t)\ < max(||x||", llyll*^) and \h{Q,y,t)\ < e~^*||y||". This is enough to show (j27|) . 

We notice that h{x,y,t) — )■ Yli=i^i^ point-wise as i — +oo. This, together with the Lebesgue 
dominated convergence yields 

, lirn^ .^■Tji.) = E = E (i.. I) = (.) 

The calculations above are valid for any choice of u, in particular, letting u \ we obtain 

(*) =^x,^/,^^ =^Xi{^^,99^ =xo(grad/,c^). 

For the sake of notational simplicity (|28p will be proved for d = 1. By the above calculations 

e^'-'+^^Tt+uRx) - X (/', if) = e>^^ {e^'TMx) - x (C, ^)) . 
Let us now treat the inner expression using (jlOp : 

e^'TMx) - X (4, (^) = E (e^* (/„(xe-^* + ou{t)G) - /„(G)) - xUG)) 

= Efe^*y^ lUy)dy-xl'JG)\ 



E(e^*y^ (C(y)-C(G))dyj + e^\ou{t) - 1)KGUG) . 
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The second term is easily upper-bounded by Ce for some C > 0. The first one can be rewritten as 
Ee'^* / / l'lXz)dzdy <^e^'\xe-^'* + {ou{t) - l)Gf sup \l'^{z)\ = {*) . 

Jg Jg zG[G,a:e-Mt+o«{t)G] 

By the discussion at the beginning of the proof we know that I" grows polynomially, i.e. there exists n 
such that 

(*) < e~^*E(x + e^\ou{t) - l)Gf ma^{\\x\\ + ||G||)" < e~^'\l + ||x|r). 

This is enough to prove (|28p . 

Now, we need an estimation of e^^Tt{lu — f){x) which takes into account u. Using the same trick as 
before it is enough to prove an estimation for e^*Tt{lu+i — h){x)- Denoting ku{x) = lu+i{x) — li{x) we 
have 

e^'TtK{x) = I e'^'iKixe-'^' + ou{t)y) - K{y)My)dy. 
Using the mean value theorem we get 

h{x, y, t) := e'^*(/„(xe-^* + ou{t)y) - lu{y)) = {x + e^\ou{t) - l)y) o grad lu{xo). 

Further following arguments used in (j30p and (j3ip and after them we conclude that in order to obtain (p9|) 
we need to estimate grad(/u+i — li). For simplicity we will provide the details for d = 1. We have 

{{lu+i - h){x)y =( [ {h{x - y) - h{x))g^{y)dy) = [ {l[{x - y) - I' {x))g^{y)dy = (*). 
Using the mean value theorem there exists a function xo{x,y) such that xo{x,y) £ [x — y,x] and 

(*) = / yli{M^,y))9u{y)dy = {**). 

JR 

The function /'/ is grows at most polynomially. Therefore for some c, n we have 

|(**)| < c / |y|(l + ||x|r + \\y\n9uiy)dy < c(n)(l + ||x|r), 

jR-i 

where c(u) is some function fulfilling the required conditions. □ 

Let us denote Yt{f) := Ft {J^Xtj) - \Xt\ (^f ,^)) , where Ft = e'^^''/^)* when Ap < 2/x; Ft = 
e-(^p/2)ti-i/2 -^vhen Ap = 2/x; Ft = e~^'^^~^^'>* when Ap > 2/i. Let us assume that we proved 

hmsupE - Yt{lu)f = c(n), (32) 

t— >+oo 

where c{u) is some function such that c{u) \ as u \ 0. Let L(lu), L{f) be the laws of the limits asserted 
in Theorem 13.21 and Theorem 13.61 

Assume that we know already that convergences in Theorem 13. 21 and Theorem 13 . 6 1 hold for any lu,u > 0. 
One can check that in either case 

m{L{lu),L{f)) ^0, as n \ 0. 
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Let us fix e > and choose u > such that c(n) < and m{L{lu) , L{f)) < e. We can choose Ti > such 
that for any t > Ti we have m{L{lu), C(Yt{lu))) < e (where C denotes the law of given random variable). 
Let T2 be such that for any t > T2 the L2 norm in (|32|) is smaller than 4e^. Now, for t > max(Ti,T2) we 
have m[L{f),Yt{f)) < 4e and therefore all of theorems mentioned hold also for /. The case of Theorem 
13.101 follows similarly, but one has to use convergences in probability and the L^-norm. 

Let us sum up this section. In the proofs below we may work with additional assumption that / is a 
smooth function and its derivatives grow at most polynomially fast if only we show (|32p . 

4.5 LLN and CLT for small branching rate 

First we prove the following fact 

Fact 4.4. Let {Xt}^yQ be the OU branching system and Xp < 2^ and let f & V . Then 

Ex (e-^^''/^)* (^Xt, /)) ^ ast^ +00. 

(e-(^^/2)* {X,j)y ^ a}, Var. [e~^'^/'^' {x„ f)) ^ a}, (33) 
. Moreover, 

supE, (^e-(^^/2)i (^(^Xt,f)-K^ (^*'-^~)))^ < +°°- 
Proof. First we note that by (I22p and Lemma [4.31 (ineq. (j27p ). 

\w'{x, t,0)\< e(^p-^)*(l + ||x|r) < e'-^p/^^^l + Hxf), 
which, by (|23p . implies the first assertion. Using ([22]) and (|23p again we calculate the second moment 



where ctj is the same as in 



{x)ds 



By (j27p in Lemma [4.31 the integrand in the last expression can be estimated as follows 



Ttf\x) + 2\p / Tt- 



(x)ds. 



Tt 



t-s 



,{Xp/2)s 



Tsh-) 



(x) < e^^^-'^>Tt-s [(1 + 



\n\21 



[X . 



Using representation (jlSp it can be checked that for any t > we have % [(1 + || • ||")^] {x) ^ (1 + 



X 



|2n\ 



The dominated Lebesgue theorem implies (|33p . We also conclude that for any t > 0, 

«;"(x,t,0) <eV(l + ||xf"). 
Similarly we investigate w'" {x,t,0). By (|22p we have 

\w"'{x,t,0)\<e^^'Tt\ff{x)+ f e^-'^'-'^Tt-s ^ (•, 0)u;'(-, 0)] {x)ds. 

Jo 
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Using the estimates obtained already and the fact that 7j [(1 + || • ||")^] (x) ^ (1 + 

\w"'{x, t, 0)1 < (1 + ||xf ")e^*'* + e^"' f e^^''^'Tt-s [l + || • f "] {x)ds < e^'^/^'^^^\l + ||xf 

Jo 

Finahy, we wiU also need the fourth moment. By ()22p and the estimates above we get 

< e-^^'Ttj\x) + e-^^^' f e^^^'-^'^Tt-s [w"{; s, + w"'{; s, 0)w'{; s, 0)] {x)ds 

Jo 

< e-^^'il + ") + / e^^'Tt-s [(1 + || • f ")] {x)ds < (1 + 

Jo 

It is now easy to get the last assertion of the fact. □ 

Now we will prove representation ([10]). We]iaveTsHi^,i2,...,ia{x) = e"(*i+*2+...+id)Miij.^ .^(-j;). There- 
fore 

I ^ fii,i2,..;id^h,i2,--;id 

\ \ii=0,i2=0,...,id=0 



° \ \ n=0,i2=0,...,id=0 

+00 +00 

= E 

+ 2Ap fiui2,-,id I 

il=0,i2=0,...,jd=0 11=0,42=0,. ..,id=0 

Now (jlOp follows easily. We are ready for 




Proof of Theorem \3.1\ (sketch). The proof of [101 Theorem 6] can be checked to hold also for the branching 
mechanism introduced in our paper. We will show now that the almost sure convergence holds also for 
bounded functions which are not compactly supported. We decompose / = f^ — f~, where f~^{x) = 
/(x)l{j(2.)>0} and f^{x) = —f{x)ls^f(^x)<o}- 1^ enough to prove the claim separately for /"*" and /~, hence 
we assume that / > 0. For n G N we consider functions hn{x) := min(max(n — x, 0), 1) and Qn : i— t- [0, 1] 
given by gn{x) := /in(||3;||). By |10| Theorem 6] we know that e~^p^ {Xt,gn) — Voo {gn,^) as., we know 
also that e~'^p^\Xt\ — )• Voo a.s. Therefore e~^p^ {Xt, (1 — gn)) — )■ V^oo (1 — dn)) a.s. Now we estimate 

{Xt, fgn) < {Xt, f) < {Xt, fgn) + ll/lloo {Xt, (1 - gn)) ■ 

Using the previous considerations and the fact that fgn has compact support we get with probability one 

(/5n,V5)Ko <liminfe-V (X*, /) < limsup (x*,/) < (/<7„, y^) Foo + H/IL ((1 " 5n), ^oo- 
t/-+oo ty+oo 

To conclude, we observe that {fgn, (f) {f, ip) and {fgn, ip) + ||/||oo ((1 - 9n),^) (/, ^) as n +oo. 
Let now f £ V and f ■= f — (/, f). Calculating similarly as in the proof of Fact 14.71 we get 

2" 



et := (e-^r>t ^Xt, /))' < e-^^'TtP{x) + ^ e'^^'Ts 



Tt-sfi-: 



{x)ds. (34) 
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Obviously for any x G we have Ttf{x) — t- 0, moreover |7t/(x)| ^ 1 + Standard considerations 

using the Lebesgue dominated convergence theorem yield that et — )• 0. Therefore e~^p^ (^Xt, f"^ — 0, 
and further e"'^^* {Xt, f) — e~''*p*|Xt| (/, if) —?-^ 0, which concludes the proof. □ 

In the proofs below we will use the term "X is asymptotically equivalent Y" to denote the situation 
that — It — 7- in probability (equivalently in law) as t — t- +oo. Now we are ready for 



Proof of Theorem 13.21 We start with the following random vector 

Zi{t) := (^e-^^'\Xtle-^^^/^^\\Xt\ - eVy^), e-^^^-Z^)* (x*,/)) . 



Let n G N to be fixed later and let us write 



where denotes the subsystem originating from the particle Xs{i). We know that with probability 

one e~^p^\Xt\ — t- Voq. Therefore e~"'^p^\Xnt\ — e~'^''^\Xt\ — t- a.s. as t — t- +oo. Let us consider the second 
term 



\Xnt\ - e^^^'Voo = \Xnt\ - e"^''* lim e'^f |X 



i,nt I 



nt+s\ = \Xnt\ - e"^-* lim e-^-(^+"*) V |X 

i=l 

\^nt\ \Xnt\ 

\Xnt\ - Hm e-'nxr'\ = ^ " 



i=l 



i=l 



where are independent copies of Voo (note that formally they depend on t, however we will suppress 
this fact in the notation). We couple each X*'* with the branching system starting from one particle 
located at 0. To this end we use the same methods as in Fact 14.11 for particles movements and retain the 
branching structure. The coupled system is denoted by X*'*. Let us write 



H 



i=l 



< 



i=l 

l^(n-l)tl 



1=1 



i=l 



By the approximation argument from Section [4.41 we may assume that / and its derivative can be upper- 
bounded by a polynomial of order A; E N. Indeed ([32|) holds by (p3|) and easy calculations. 



18 



By FactOand the mean value theorem we have |/(X[;*_i)j(j))-/(^[;*_i)i(j))|l||Xt{j)||<t ^ te 
II (n-i)t(-^)ll)''' Using the conditional expectation with respect to Xj, (|23p and 7t(l + 1| ■ ^ (l+||x||'^) 

we get 

i=l 

+ e-("V2)tE,.^l|l^^(^)ll>,(l + ||X,(i)||)^e("-i)V 

<g(nAp/2)t^g-(n-l)M ^ g((n-2)Ap/2)i-7^((-^ ^ || . \\)''^^^^^)(x) . 

There exists uq > such that for any n > no we have nXp < 2^{n — 1). Using the Schwarz inequality, the 
second term can be estimated by a/7((1 + || • ||*^)^(rc) Y^7tl||.||>((x). The Ornstein-Uhlenbeck process has 
Gaussian marginals with bounded mean and variance therefore 7tl||.||>((a;) < e"'^*^ for a certain c > 0. We 
may conclude that E^^// — ?• as t — t- +oo. We notice that X*'*'s are i.i.d. branching particle systems; indeed 
the only dependence among X*'* is by the initial condition. Let us denote Zl := e~^^'^~^^^^^'^^^ / 
and zl := EqZ'I (we will write Eq to underline the fact that X*'* starts from 0). 
By (p3l) and Lemma 14.31 one checks easily that for n large enough 

E^e-(Ap/2)t ^ 1^*1 ^ e(Ap/2)t|E^^l| < g(nA,/2)tg-(n-lW ^ ^ ^ 

Using the facts above we conclude that there exists n such that Zi{nt) is asymptotically equivalent to 

/ \Xnt\ \Xt\ 

Z2{t) = e"V|Xt|,e-("^''/2)t ^ (1 - V^) ,e-(^''/2)i ^(^i _ 

\ 1=1 i=l 



Let us now denote 



(\Xnt\ \Xt\ \ 

e-V|x,|,|X„,|-'/'E (l-^~)'l^*r'^'E(^*-^*) ' 
i=l 1=1 ) 

which we will consider conditionally on the event {\Xt\ 7^ 0} (here and below we adopt the convention 
that a/0 = 0, Z^ is well defined then). The corresponding expected value is denoted by E^^ |Xi|^o- Let us 
denote the characteristic function of Z^ 

Xi{0i,02,es;t) := 

( \Xnt\ \Xt\ 1 

E.,|x,|^o exp I ie,e-^-'\Xt\ + iO^lX^tr'^' ^ (l - + i^sl^r'^' ^(^t " ( ■ 

[ i=l i=l J 

Conditioning on Xnt and using the Markov property we check that variables 1 — are i.i.d, moreover 
they are independent of the system before time nt. We denote their common characteristic function by h. 
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We have 



ieie-^^'\Xt\ + ies\Xt\-'/^ Y.^z\ - zi) \ h (e^/VW^ 



1 = 1 



By Fact 12. ll and the central hmit theorem we write h{02/ ^/n)"^ — )• e ^l/'y'^'^v) ^ where cjy 
we will work with 



V ~ 2p-l 



. Further 



X2{ei,92,e-y,t) :-- 



\xt\ 



The limit of xi is the same as the one of X2 providing that any of them exists. Indeed 

\x2{9i,02,e3;t)-xi{0i,e2,e3;t)\ 

h[e2/^Je(^-^)^p'\Xt\) -h(92/^Ax^^ 



„(n-l)Apt|jj'^ 



The sequence of events {\Xt\ 7^ 0} decreases to Ext^ and ¥ (Ext^) > 0. We have \h\ < 1 hence the first 
summand converges to 0. By Fact l2.1l on Exf^ we have 



\X„t\ 



|jSfj|g(ri-l)Apt 



la.s. therefore the second summand 



converges to as well. Now we recall that are i.i.d and independent of Xt- Conditioning with respect 
to \Xt\ we get 



X2{9i,e2,e3;t) = 



g{n-l)Apt|j(^^ 



(35) 



where h denotes the characteristic function of Z^ — z} . We will now investigate the second term of the 
product above. To this end we fix two sequences {tm} and {am} such that tm — +00, — ?• +00, as m — >■ 
+00 and consider 

1=1 

where Zl^ are i.i.d. copies as described above and := flm Varo(.^t^). From Fact 14.41 we easily see that 
is of the same order as a^- We will now verify the Lindeberg condition. Let us calculate 



-1 (^m 



i=l 



^1 \2i 



By by the Schwarz inequality, the Chebyshev inequality and Fact 14.41 we get 
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as m — 7- +00. The CLT holds therefore Sm. -^'^ -^(0, 1). Obviously ^ ^ (^"j- We conclude that 

hz,^^..^ {es/V^r- ^ e-^3/(2.^). (36) 

Now we go back the X2 

X2(^l,^2,^3;*) = ^x,\Xt\^oi- ■ ■)'^Ext +'^x,\Xt\i^o{- ■ ■)'^Ext''- 
Arguing as before we get that the first summand converges to 0. The second one is 

IEx',|Xt|^o(- ■ ■)'^Ext'^ = '^x,\Xt\^o{Ext'')'Kx^Ext''{- ■ •)) 

where '^x,Ext'' denotes the conditional probability with respect to Ext^. We have ¥^Xt\^oi^^'l^'^) ~^ 1- 
We proved (j36p for arbitrary sequences therefore the second term in ()35p converges (on the set of non- 
extinction) as follows 

(03/1^*1'/')"''' ^ e-'"/'-'''f\ a.s., as t ^ +00. 

We know that e~'^''*'\Xt\ — s- W a.s. (recall that W denotes Voo conditioned on Exf^). By the Lebesgue 
dominated convergence theorem we get 

X2(^i, ^2, ^3; t) e^^2/(2a^)e"^3/(2<xp]£^ i^iQ^w} , ast^ +00. 

We know that {\Xt\ 7^ 0} is a decreasing sequence of events and Hi^iil^tl 7^ 0} = Ext^. Moreover 
P {Exf^) > 0, therefore by Lemma [4.21 the limit of Z^^t) is the same if considered conditioned with respect 
to {\Xt\ ^ 0} or conditioned with respect to Exf^. In this way we have proved that conditionally on Exf^ 

where (VF,Gi,G2) is the same as in the statement of Theorem 13.21 Using some obvious transformations 
we get 

Z2{t) {w, ^GuVWG2) and Zi{t) -^'^ {w, ^GuVWG2) , 
which easily implies the convergence asserted in the theorem. 

4.6 CLT for large branching rate 

We start with 

Proof of Fact [3731 The fact that H is a martingale with respect to the filtration of the OU branching 
system follows by its branching property and easy calculations using (j22p . (|23p and (jl9p . Consider now 
f{x) = \\x\\. Using ([22]) and (f23|) we can get 

Jo L J 

< e(-^f+2A')*(l + \\xf) + e(-^f+2M)t f e^^'Tt-s [e''^'"' f] {x)ds 

Jo 

< (1 + \\xf) (1 + e^-^^'+^/^^^^-^^dsj < C{1 + \\xf). 
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To check that V and H are orthogonal one can use (l23l) together with the polarisation formula. We leave 
this task to the reader, just mentioning that this result is somehow expected as l,x are orthogonal with 
respect to densities gt defining (jlSp . □ 

Proof of Fact COi The proof follows easily by Fact 14.1 1 Indeed it is enough to notice that, if we couple X 
and X according to the procedures described above, we may write 

Ht-Ht = e(-^+'^)*xe-^*|Xt| = e-^^x\Xt\, 

where H,H are defined according to (|14p . □ 

It will be useful to know that the range of the OU branching system grows at most linearly. This is 
a well-known fact fact for the branching Brownian motion. It is somehow obvious that it also holds for 
the OU branching motion as the Ornstein-Uhlenbeck process is "better concentrated" than the Brownian 
motion. However, we were not able to find a proof in the literature hence we provide one. 

Lemma 4.5. Let X be the OU branching system starting from x. There exists a constant C and a random 
variable T such that with probability one 

Vt>TV,6{i,2,...,|x,|}^t(i) G i?(xe-^*,Ct), 
where B{0,r) denotes a ball of radius r centred at 0. 

Proof. By the usual coupling argument it is enough to prove the fact for x = 0. Let us denote by An an 
event that X up to time n is contained in B{0,C{n — 1)) for some constant C. By the Borel-Cantelli 
lemma to show our claim it is enough to prove X^„>]^ IP (A'^) < +oo. By the result from |16| . the coupling 
construction and Gaussian concentration inequality we know that the supremum of the Ornstein-Uhlenbeck 



process on the interval [0, n] can be stochastically dominated by Ci(ylog(n) -|- where Ci > is a 
certain constant and G is a standard normal random variable. Let us fix any 7 > Ap and estimate 



«) < P K n {\Xn\ < e^"}) +P(|X„| > e^") 

= P( 3i6{i,2,...,n} sup \\Xs{i)\\ >C(n-l)n{|Xj <eT"} ) +P(|X„| 
V se[o,n] J 

+ |G|) > C{n - 1)) e^" + P (|X„| > e^") , 



< 



where Xs{i) denotes the position of the particle i or its ancestor at time s E [0,n]. The proof is concluded 
by using the estimations of the Gaussian tail in the first term and Fact 12.11 together with the Chebyshev 
inequality for handling the second one. □ 

We will also need 
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Fact 4.6. Let {Xt}^^Q be the OU branching system starting from x and Xp > 2fi. Moreover, let f G V . 
Then 

supE, (^e-(^^-^)* (Xt, /))' < +00. 
Moreover there exists c : M+ 1— )■ M+ such that c{u) — )• as n \ and 

limsupE, (^e-(^-"^)* ((Xt, /) - (X*, < c{u), (37) 

where 1^ is given by (|26p . 

Proo/. We recall ([23]) and By Lemma O one sees that w'{x,t,{)) < (1 + ||3;||")e('^''-'')*. Now we 

check that 

ft 



(x)ds 



w"{x, t, 0) < e^''*7I/2(x) + ^ e^''^*-^)?^-, ((1 + || • \\^)e^^p~^^ 

< e^^'TtP{x) + e^-^' f e^^-^-'^^^'Tt-s [(1 + || • f")] (2:)ds < e^^^-^-^''^\l + 
Inequality (|37|) follows by (|29p and the above calculations. Indeed, one checks that 

limsupe-2(-^f-^)*t(;"(x,t,0) < c(n)(l + ||x||"), 
for some function c(n) such that c{u) — )■ as n \ 0. □ 

Proof of Theorem 13.101 Our first aim will be to prove the convergence of the spatial fluctuations. To 
this end we denote 

Yr{t) := e-(^--^)* {{Xt, f) - \Xt\ (/, (^)) = e~^^^~'^^' {Xt, f) , 
where f{x) = f{x) — (/, 92). We have 

\Xt+s\ \Xt\ \x'/\ 

yi(t + s) = e-(^''-'^)(*+^) ^ =e-(^*'-^)* J^^e-^^*--^)^ ^ f{Xi'\j)) 

i=l i=l J=l 

\Xt\ /\xl-'\ 

i=l \ j=l 

where denotes the subsystem originating from the particle Xt{i) and {Xl'^}s is a coupled version 

of X*'* starting from 0. We write 

\Xt\ \xi''\ \Xt\ 

Y,{t + s) = e-("--^)* (/W*(J)) - + e-(^--'^)(*+^) (x^'*, /) . (38) 

i=l i=l 2=1 

Let us first deal with the second term 



i=l 
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where fs{x) := f{x) - 7^/(0) and Ig := 7^/(0) - {f,(p). By Lemma [43] we know that < e"^'''* and 
therefore 

From now on, whenever we prove convergence, we will assume that s = 2t. With this convention the 
above expression converges to 0. We denote the first summand of 12 by 13. As we explained in the proof 
of Theorem 13.21 the systems X*'* are i.i.d. Moreover their only connection with Xt is via the number of 
particles \Xt\. Below, we will use the conditional expectation given Xt which will be denoted as To 
ease the notation we will write simply EqX*'* instead of Ej^^X*'*. Using the conditional expectation and 
the particular choice of fs we calculate 

\Xt\ \Xt\ 

(Y.it + s)f = e-2(^p-^)(*+^)E, Yl E ((^^'*' f^) {^s\ fs)) 

i=i j=i 

= e-2(^''-'^)(*+%.f;Eo((x:'*, 

By Fact 1321 one can prove that Eq (^(^Xi'\ /^^^^ < e^^^f"^)^ therefore 

E, (y3(t + s))2<e(-Ap+2M)t^0_ 

In this way we proved that Y2{s + t) — t- in probability. By arguments in Section [4.41 we may assume that 
/ has the second derivative which is bounded by a polynomial. To this end we check that (j32p follows by 

dSZD. 

Now we decompose the first term of (|38p according to the formula 

Y5{t + s) + Ye{t + s) := 

\Xt\ \xi'*\ \Xt\ \xi'\ 

g-(A,-M)(i+.) ^ ^ i?(i,t,s) + e-(^''-^)(*+^)^ ^ gradf{X^\j))o{X^\j)-X^\j)), (39) 

i=l j=l i=l j = l 

where E{i,t,s) := - f{Xi'\j)) - gr ad f {Xi'\j)) o {xi'\j) - Xi'\j)). We now notice that by 

the coupling properties we have 

Xl^\j)-Xl'\j) = XS)e-'^^. (40) 

We have \E{i,t,s)\ < R''{\\Xt{i)\\e->''f if only xi'\j) , xi'\j) belong to the bah of radius R. Let us now 
denote the event that this condition holds for all i,j with R = C{t + s) + ||x||e~^* by At- By Fact I2.ll we 
thus prove 

E,|y5(t + s)\1a, < e^(^+*)((t + s) + e-^ixID^^e-^/^^ ^ 0. 
By Lemma [4.51 we have P (At) — )• 1, hence we proved that Y5(t + s) — )• in probability. Using again (1^0]) 
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we may rewrite Ig as 

Yeit + s) = e-^^''-^)* ^ Mi) o {Zi -Zs)+Zsol e-(^^-^)* ^ Xt(i) 

i=l \ i=l 

( \xt\ \ I \Xt\ 

+ (grad/,99)o ( e-(^--^)*^Xt(i)(e-^-^|X:'*| -1) 1 +(grad/,^)o I e-(^^'-A')* ^ 

= : Yj{t + s) + Y^{t + s)+ Y^{t + s) + YiQ{t + s), 

where h{x) = grad/(x) — (grad /,(/?), := e"'^^** and := EqZ] (this does not depend on i 

since Z] are i.i.d.). Simple calculations using <\22} and Lemma [4.31 (second estimation in (j27p ) reveal that 
ll-^sll ^ e~'^^^ . Moreover by Fact 14.61 the covariance matrix Covo(Z]) is bounded by some C (in a sense 
that each entry is bounded). Using conditioning with respect to \Xt\ we have (we skip the index ) 



W7{t + sf = e-2(^^-A')*E, ^*(^)^ Covo {Zl)Xt{i) < e-2(^^-^)*E, ^ \\Xt{i)f ^ 0, 

1=1 1=1 

where ^ denotes the transposition and the convergence holds by assumption that Ap > 2/i. The conver- 
gence: K.j,{Yg{t + s))2 ^ follows in a very similar fashion and is left to the reader. Now by Fact 13.81 
one easily checks that l8(^ + s) — s- a.s and Yiq converges to the same limit as in the thesis of Theorem 
I3.1UI A revision of the steps above reveals that the limit (in probability) of Yi is the same as the limit of 
YiQ. This in fact is also enough to conclude the second convergence in Theorem 13.101 To end the proof we 
notice that to obtain the joint convergence (jlSp one may use the same methods as in the proof of Theorem 



4.7 CLT for critical branching rate 

Fact 4.7. Let {Xt}^^^ be the OU branching system and Xp = 2^. Moreover, let f £ V , then there exists 
a constant C such that 

(e-(^''/2)i^-i/2 ^;s^^^ ^0 ^ ^ 



E 



[e~^\nH-y^ {X„ ^ aj, Var. (^e~('./^^H^y^ {x„ /)) ^ a% 



(41) 



where cjj is the same as in (fT2|) . Moreover, sup^^E^^^ ^/^ i^Xt, /yj < +oo for any 6 > 0. 

Proof. The first convergence follows easily by (j27p in Lemma 14.31 Using (I22p and (|23p again we calculate 
the second moment 



E, (e-(^-/2)t^-i/2 ^^^^ j'^y ^ t-^Ttf\x) + 2Xpe-^-H-^ e^^^'-'^Tts {e^^'TsR 



t"^TtP{x) + 2\pt~^ frt- 
Jo 



Tsfi- 



{x)ds 
{x)ds. (42) 
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We recall that Ap/2 = fi. Using Lemma [4.31 (equation (I28p ) and elementary considerations, we obtain 
that the limit of the above expression is the same as the one of 




One easily checks that Tt[xiXj]{x) — )■ hence the second term disappears in the limit. We also have 
7t[x?](x) — )• the whole expression converges to aj given by (I12p . 

Obviously the limit of variances is the same. We also conclude that for any t > 

w"{x,t,0) < (1 + ||xf")e^f*t. 

Similarly we investigate w'" {x,t,0). By (|22p we have 

\w"'{x, t,0)\< e^^'Tt\ff{x) + f e^-(*-^)7I-. [w"{; s, Q)w'{; s, 0)] {x)ds. 

Jo 

Using the fact that by (f27|) |'u;'(rE, t, 0)| < e^'^p~^^^{l + HxH"), together with the above estimate on w" and 
the fact that % [(1 + || • (x) < (1 + ||xp"), we get 

\w"'{x, t, 0)1 < (1 + "■)e^''* + e^^* f e'^^^/^> sTt^s [l + II • f "] {x)ds < e((3/2)Ap)ti(i + ||^||3n)_ 

Finally, we will also need the fourth moment. By (|22p and the estimates above we get 

+ e"2Vr2 r e^p'^'-'^Tt-s [w"i-, s, Of + w"'{-, s, 0)w'{-, s, 0)] (x)ds 
Jo 

< e-^^H-^l + ||xf ") + e-^^H-^ [ e^^'s^Tt-s [(1 + || • f ")] {x)ds < (1 + ||xf 

Jo 

It is now easy to check that for t > 6, 

(e-(^^/2)*rV2 ((x,,/) - E,. {Xt,f))y < (1 + llxf"). 

□ 

We have yet to prove (fT3|) . By ([7|) we have -£^(p{x) = —2{n/a'^)xiip{x). Therefore, 

(^'^) = (^'S) = (/'-2W^')^i'/'> = -2W^')/i,o,...,o. 

Other coordinates can be treated in the same way, giving ()13p . 
Now we are ready for 
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Proof of Theorem 13. 6L In this proof we will us both ideas of the proof of Theorem 13.21 and Theorem 
13.101 We start with the following random vector 

Let n G N be fixed later and let us write 

where denotes the subsystem originating from the particle Xs{i). Analogously in the proof of 

Theorem 13.21 the second term is equal to e-("-*^p/2)t ^|^!|^'' (l — V^), where are independent copies of 
Voo arising from i-th particle. Next, we couple each X*'* with the branching system starting from one 
particle located at 0. To this end we use the same methods as in Fact 14.11 for particles movements and 
retain the branching structure. The coupled system is denoted by X*'*. We write 

/\Xt\ \Xt\ 
^1=1 i=l 

By arguments in Section 14.41 we may assume that / has the second derivative which is bounded by a 
polynomial. To this end we check that (|32p follows by (j4ip (and some elementary calculations taking into 
account the form of the limiting variance) . This expression is harder to analyse compared to the case of the 
small branching rate. Now we will use the methods and notation of the proof of Theorem 13.101 Recalling 
that Xp = 2/i one notices that HJ^ is the same expression as (nt)^^/"^ (^(^ + {n — l)t) + YQ{t + (n — l)t)) 
when one puts s = {n — l)t (see (p9l) ). Is can be handled in the same way, if only n > 2. Next, using the 
coupling properties we decompose Ig in the following way 

\Xt\ \Xt\ 
i=l i=l 

+ ksO ||(nt)-V2e-(A./2)*^^^(,)j ^ (43) 

where this time Kl := e""^*"^ ^grad/, X]'*^ and ks := Eg-fC] (this is independent of i as Kl are i.i.d). 
Eoll-f^s ~ S'l's uniformly bounded in s which follows in the same way as estimation of et in (fM|) . Using 
conditioning in a similar manner as in the estimation of I7 in the proof of Theorem 13.101 we obtain 

\Xt\ \ ^ \Xt\ 



I (nt)-i/2e-(Ap/2)t^Xi(i)o(i('^-fc,) =(nt)-ie-^-*E,J]Xi(i)^Covo(i^':)Xi(i) 

i=l / 1=1 

\Xt\ 

< (nt)-ie-VE, ^ \\Xt{i)f < {ntr\l + \\x\n ^ 0, 



1=1 
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where again we used (l22l) and (I23p . Further, we notice also that {x,Lp) = therefore by Fact 14.71 and 
the fact that kg is bounded we know that the second moment of the second term of (|43p is bounded by 
c/{2^/n), where c > is a certain constant. Thus by ()25p . we may conclude that for any n there exists t„ 
such that for any t >tn, 

ll^rib < 2c/ V^. (44) 
We recall that X*'*'s are i.i.d. branching particle systems. Let us denote 

1 /2 

and zt := EoZi By Lemma |43] one checks that \zt\ = e(("-i)^p/2)t(j^^)-i/2|^j^Q^)| < g{(n~i)(Aj,/2-2/.))t^ 
Therefore for n > 2 we have 

Ee-{Ap/2)t ^ ^ e{Ap/2)t|jE^^l| < g(A,/2)tg{(n-l){A,/2-2M))t ^ ^ ^ 

1=1 

Using the facts above we conclude that for any fixed n > 2, the expression Zi is asymptotically equivalent 
to 

i=l 1=1 

This expression differs from the analogous in the proof of Theorem 13.21 onlv by additional term Let 
us skip it for a moment and write 

/ l^ntl \Xt\ 

y i=l i=l 

Returning to the flow of the proof of Theorem 13.21 one can prove that 

1/2 



conditionally on Ext^. In order to achieve it one simply uses Fact 14.71 instead of Fact 14. 4[ The additional 
term (^^^)^^^ originates from the same term in (j45p . 

To finish the proof let us denote the law of the triple in the last expression by £„ and the one of 
{W,VWGi,VWG2) by Coo- For any e > 0, there exists n such that m(£„,£oo) ^ e and cj^fn < 
e, where c is the constant in (|44p . Now we choose T large enough to have m{C{Z'!^[t)),Cn) < e, 
m{C{Z^{t)),Z^'{t))) < e and m{C{Z^{t)), C{Zi{nt))) < e for any t > T. By applying the triangle in- 
equality we get m{C, C{Zi{nt))) < 4e for any t >T. The proof is concluded by obvious transformations. 
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Appendix 



Proof of Fact [^T/l Using (f20|) with f{x) = 1 we get (we drop argument x) 

^w{t, 6) = {\p)w{t, ef - Xw{t, 6) + A(l - p), w{0, 9) = e"^. 

The solution of this equation is (see e.g. jj] Section III. 5]) 

A(l - p){e-^ - 1) - e-^^\\pe~^ - A(l - p)) 
^~ Ap(e-^-l)-e-V(Ape-^-A(l-p)) ' 

Now we want to investigate convergence of e~'^^^\Xt\. Its Laplace transform is 

^ A(l - p){e-'^-'^' - 1) - e-V(Ape-^-"^''' - A(l - p)) 
^ Xpie-o^'^"' - 1) - e-^vt[\pe-de-^^' - A(l - p)) ' 

Using the first order Taylor expansion and dropping terms of lower order we get 

A(l - p) {-9e~^pt + o{ee~^p^)) - e~^p^ {-XpOe'^^^ + o(e-^f*) + Ap) 



L{t,d) 



-A(l -p)(9e-V _ \pe-^vt 



—Xp9e — XpC ^p^ 

Therefore 



Lit, 9) \ = ^ -—^ =: L(9) as t +oo. 

^ ' XpO + Xp +p{2 + 0)-l ^ ^ 



Taking 9 — )• +oo it is easy to check that ¥ (Voo =0) = Pe (we recall that pe = -^), therefore Voo > on 



the set of non-extinction Ext^. Let us now calculate the law of on the set of non-extinction 
L{e) = Ee-^^- = Ke-^^^lExt + Ee-^^^lExt^ = + (1 - p,)E (e-^^-|Sxf ) . 

Therefore 



E ( e-^^°°\Exf 



Xp {2p - 1) 



Xp9 + Xp p9 + {2p - 1) ■ 



Further 



4 _ d'^w{t, 0) _ e^^p (-1 + 2 (-4 + 7e*^f ) p + (8 + 16e*^f - SGe^^^p) p^ + 8e*^p (-2 + Se^^^f) p^) 
d9^ (-l + 2p)3 ■ 

Now the second part of the fact follows. In order to prove that all moments are finite we notice that 
derivatives of ()46p are of the form 



d0" ' (Ap(e"^ - 1) - e-V(Ape-^ - A(l - p)))2" ' 

where l{t,9) is a certain expression. Obviously the denominator is finite for 9 = hence the proof is 
concluded by the properties of the Laplace transform (e.g. |12| Chapter XIII. 2]). □ 
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